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Abstract 

Linear global modes, which are time-harmonic solutions with vanishing boundary 
conditions, are analysed in the context of the complex Ginzburg-Landau equation with 
slowly varying coefficients in doubly inhnite domains. The most unstable modes are 
shown to be characterized by the geometry of their Stokes line network: they are found 
to generically correspond to a configuration with two turning points issued from opposite 
sides of the real axis which are either merged or connected by a common Stokes line. A 
region of local absolute instability is also demonstrated to be a necessary condition for 
the existence of unstable global modes. 

1 Introduction 

Wakes, jets, counterflow mixing layers are well-known to exhibit under certain flow 
conditions self-sustained oscillations. It is well-established that these vortical structures 
result from a Hopf bifurcation and are due to a global destabilization of the underlying 
non-parallel basic flow (Mathis et al. 1984, Monkewitz et al. 1990, Strykowski Sz Niccum 
1991). Furthermore, they seem to be connected to the appearance of a sufficiently large 
region of local absolute instability (Monkewitz 1990). Most theoretical investigations 
have so far used a WKBJ approach which assumes the basic flow to be weakly non¬ 
parallel in the streamwise direction. As explained in Huerre & Monkewitz (1990), 
it permits to define local stability properties and to link the global spatio-temporal 
behavior of the perturbation in the streamwise direction to the local dispersion relation. 

The complex Ginzburg-Landau equation with slowly varying coefficients appears to 
be the simplest model providing a scenario of transition consistent with experiments 
(Huerre & Monkewitz 1990). Global modes, which are time-harmonic solutions satisfy¬ 
ing vanishing boundary conditions have been analysed, either with linear or quadratic 
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coefficients (Chomaz et al. 1988), or under specific assumptions on the number of turn¬ 
ing points involved in the frequency selection criterion (Chomaz et al. 1991, Le Dizes 
et al. 1994). 

The main objective of the present article is to obtain a characterization of the most 
unstable global modes without any assumption on the coefficients. An upper bound for 
the growth rate Am[a;] of the most unstable global mode is obtained in section 2. In 
section 3, the evolution with respect to the frequency of separation curves -or Stokes 
lines- is analysed in the complex plane in order to identify most unstable global modes 
conhgurations. Characteristic properties of global modes and expressions satished by 
their frequency are given in section 4. The last section recalls the main results and 
briefly discusses their extension to the weakly nonlinear regime. 


2 Necessary conditions for the existence of a global 
mode 


In the present framework, a global mode of global frequency a; is a time-harmonic 
function of the form 


^{x,t;uj,e) = ^jJ{X;u!,e)e , 


( 1 ) 


where 'll! satishes the Ginzburg-Landau equation 


1 


(Ai) fcg ( A^) 


A 

dx 


-u:uk{X)kl{X)+u:o{X)) 


ip{X;uj, e:) = 0 


( 2 ) 

and vanishing boundary conditions at -|-oo and —oo. The coefficients a;o(A), uJkk{X) 
and ko{X) are assumed to be complex analytic functions of the slow variable X = ex 
where e is a small positive parameter, and uikk{X) is nonzero in the entire complex 
plane. 

The local dispersion relation associated to has the following form : 


a; = 0 . 0 (A) 


( 3 ) 


Causality implies that there exists a local maximum growth rate for the perturbation 
of the basic flow, which means that the quantity uii^rnax{X) = meoLkeRiXiik, X)) is 
everywhere dehned as well as its maximum over all A real denoted This condition 

guarantees also the existence of and uJi^max{oo) respectively dehned as being the 
maximum of a;o,i(A) over all A real and the larger of the two limiting values taken 
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by uji^rnaxi^) at X = +00 and X = —oo. Note that the wavenumber k = k^i^X) 

du) 

has a vanishing group velocity TTr(^o) = 0. According to Bers (1983), k^iX) and 

ok 

u}o{X) = uj{ko,X) are indeed the local absolute wavenumber and the local absolute 
frequency, respectively. The sign of a;o,i(X) characterises the absolute/convective nature 
of the local instability. If a;o,i(X) > 0, the location X is said to be (locally) absolutely 
unstable. It is convectively unstable or stable, otherwise. 

It is well-known that there exist two formal solutions of (^ which admit the following 
WKBJ approximations as e ^ 0 (Bender & Orszag 1978, Wasow 1985) 

^ A^{X]uj)eel'' , (4) 


where the two local wavenumbers k^ and k are dehned as 


k^{X-,uj) 


ko{X) ± 




2 ^ - a;o(X) 
^kk{X) 


( 5 ) 


The square root can be arbitrarily dehned, but for convenience, we shall use in this 
section the value of the square root of positive imaginary part with a branch cut on the 
positive real axis. 

The solutions and are not asymptotically valid in a full complex neighborhood 
of turning points dehned]^ by 


k+{X-uj) = k-{X;uj) 


( 6 ) 


However, they are known to be uniformly valid approximations of solutions of (|2|) 
in any (bounded) domain where any two points can be connected by a path 
along which ^m[J^{k~^{S-,uj) — k~{S-,uj))dS] is a differentiable and nondecreasing 
function (Wasow 1985, Fedoriuk 1987). 

In the present context, WKBJ approximations are useful only if boundary conditions 
can be applied to them, and therefore only if they are uniformly valid from arbitrarily 
large negative X to arbitrarily large positive X. Their use then requires a precise analysis 
in the complex X plane of the separation curves ^>m[/^(fc+(S'; a;) — k~{S-,uj))dS] = Cst 
in order to identify the region of uniform validity, as well as a study of the signs of 
and near inhnity along the real axis to determine the behavior of each 

WKBJ approximation for large X. 

According to the above dehnitions, the signs of the imaginary part of both branches 
k'^ and k~ are easily determined. For all uj such that uji > and all real X there is 

^Note that these points are the branch points of the square root and satisfy also oJo{X) = lu. 
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no real k solution of the local dispersion relation (^, and when uji +oo, the branches 
and k~ are opposite. It follows that for all oji > uj'^^ax all real X, Q'm[/c’''] and 
are differentiable, of constant and opposite sign as illustrated on Figure |5(a). 
When uJi becomes inferior to oJ^maxi and '^m[k~] eventually change sign. But, 





Figure 1: Typical sketches of the imaginary part of the spatial branches k~^{X]U!) and 
k~ {X ; a;) versus the real location Xr at a hxed complex value of uj with the square root 
dehnition prescribed in the text, (a): ojY'max < (b): ‘^i,max(^) < a;* < oj^^^x'i (c): 

max[uJi^rnax{oo),UjJ^^^^] < UJi < UJ^Zx' 


by dehnition, as long as a;* > uji^maxioo), these changes do not occur near oo [Figure 
|l](b)]. Behaviors for large X of both WKBJ approximations are then invariant for all 
> ^i,max{oo) 'ip^ is exponentially small for large positive real X and exponentially 
large for large negative real X and the opposite holds for 'ip~. 

When uji > 0 ;^“, the spatial branches k'^ and k~ satisfy an additional property : 
according to the dehnition of the square root prescribed above, both branches are ana¬ 
lytic complex functions for all real X and satisfy a;) — k~{X]ui)] > 0 [Figure 

|^(c)]. The function ^m[J^{k~^{S-,u}) — k~{S-,uj))dS] is then a diherentiable and nonde¬ 
creasing function on the real X axis, so that and "0“ are uniformly valid asymptotic 
solutions in any interval of the real axis. 

The behaviors of these two independent solutions for large X hnally imply that no 
solution vanishing on both sides at 00 can be constructed Q. In other words, all global 
modes have a frequency satisfying 


UJi <max[uJi^rnax{00),UJ^Z] ' (V 

In physical terms, this in particular means that a medium which is stable at inhnity 

^The uniform validity near infinity of WKBJ approximations is not essential for the proof, although 
it can be achieved under certain conditions (Fedoriuk 1987). 
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(|^i,ma*(oo) < 0) is globally stable if there is no region of local absolute instability, i.e. 
no location where uJo^i{X) > 0. 

Note that this result is stated in Chomaz et al. (1991) with an incomplete proof and 
proved in Le Dizes et al. (1994) by another method in the case of two-turning-point 
conhgurations. 

3 Stokes line network deformations to a global mode 
configuration 

In this section, we assume that uJi^max{oo) < and show that the analysis of the 

separation curves 9'm[/'^(/c’''(S'; ci;) — k~{S-,u}))dS] = Cst in the complex plane permits 
to identify the most unstable global mode conhgurations of frequency u such that uii > 
(oo). 



Figure 2: Example of separation curve network for a frequency u such that Ui > 

If oji > the geometry close to the real axis of the separation curves is partic¬ 

ularly simple : with the dehnition of the square root prescribed above, the real X axis 
crosses each line once and only once, and branch cuts and turning points are far from 
the real axis. The domain of uniform validity of WKBJ approximations can be extended 
in the complex X plane by following the simple procedure that consists in deforming its 
frontier without cutting any branch cut, nor twice separation curves of the same label. 
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The largest domain that can be obtained is limited by the Stokes line^ that do not 
cross the real axi^ as illustrated on Figure ^ 

As long as uii > u}i^rnax{oo), none of the WKBJ approximations satishes simulta¬ 
neously both boundary conditions, as obtained above (see Figure |l|(b)). Furthermore, 
for these values, the geometry of the separation curves remains unchanged near inhnity 
because no turning points goes to inhnity and none of the separation curves become 
asymptotic to the real axis near oo. However, when uji is smaller than some turn¬ 

ing points may have crossed the real axis and some separation curves become twisted 
but the domain of uniform validity of WKBJ approximations may still connect J-cxd and 
—oo on the real axis, which guarantees the non-existence of global modes. This fails as 
soon as two turning points or two Stokes lines which were on opposite sides 
of the real axis for large oji collide. 



(a) (b) 


Figure 3: Generic pinching conhgurations ; (a): Pinching between two turning points, 
(b): Pinching between two Stokes lines. 

When such a Stokes line conhguration occurs, the domain of uniform validity becomes 
“pinched” between two turning points [Figure ^(a)] or two Stokes lines [Figure |5(b)] 
and both inhnities on the real axis are no longer in the same domain. The function 
Am[/^(/c+(S';O') — k~{S]oj))dS] becomes stationnary at the double turning point Xg in 
case (a), and on the Stokes line connecting the single turning points Xi and X 2 in case 
(b). The uniform validity of the WKBJ approximation from large negative X to large 
positive X is then no longer guaranteed and the existence of a global mode cannot be 
excluded by the former argument. 

^i.e. separation curves issued from a turning point 

redefinition of the square root may be needed to maintain branch cuts in the grey sector limited 
by Stokes lines of Figure ||. 
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One can indeed demonstrate that the frequency that corresponds to one of these 
Stokes line networks is actually the leading-order approximation to a discrete number 
of global frequencies. The next section gives a outline of the proof. 

4 Most unstable global modes characteristics 

The simplest proof is based on the local analysis of turning point regions. Since WKBJ 
approximations are still uniformly valid on each side of the region where pinching has 
occurred, the existence of a global mode is equivalent to the correct matching of the 
WKBJ approximations prescribed by the boundary conditions at inhnity^ across either 
the double turning point Xg [Figure 0(a)], or the two single turning points Xi and X 2 
connected by a Stokes line [Figure 0(b)]. 

This kind of problem is treated in several textbooks for real conhgurations (see for 
instance Bender & Orszag, 1981, chap. 10). The same analysis can easily be applied in 
the complex plane and leads to the following results. 



(a) (b) 


Figure 4: Global modes approximations in the complex X plane ; (a) : Double-turning- 
point global mode, (b) : Two-single-turning-points global mode. 

For the double-turning-point conhguration, matching is possible if and only if the 
local approximation near the double-turning-point region is a Hermite function [see 
Figure 0(a)]. This gives a discrete number of global modes whose leading-order frequency 
®i. e. both WKBJ approximations subdominant at -l-oo and — 00 , respectively. 
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is the local absolute frequency taken at the double turning point Xg : 


= a;o(X,) . 


( 8 ) 


For the two-single-turning-point conhguration, one has to proceed to a double match¬ 
ing because there are three large separated regions [see Figure §(b)]. Along the Stokes 
line connecting both turning points, any solution admits a uniform decomposition into 
both families of WKBJ waves. But, that approximation breaks down near each turn¬ 
ing point where a local approximation is built in terms of Airy functions. These local 
approximations can themselves be matched to each subdominant WKBJ approxima¬ 
tion if and only if the following integral relation is satisfied by the leading order global 
frequencies 


rX2 


Tie Jxi 


\ 


2 ^'°’ 


1 

«+2 


(9) 


^kk (^) 

This relation, where n is an integer of order 1/e, constitutes the leading-order expression 
satisfied by a discrete number of global modes. When n is of order unity, both turning 
points are distant of e and one recovers the double-turning-point case. The same results 
can also be obtained using a method based on uniform approximations as shown by Le 
Dizes et al. (1994). 


5 Discussion 

We have shown that a global mode for the Ginzburg-Landau equation has always a 
growth rate cUj smaller than max[a;j_max(c)o),Furthermore, the most unstable 
global modes of growth rate larger than uJi^rnax{oo) have been demonstrated to be de¬ 
scribed by specihc Stokes line conhgurations. They are generically[] associated to the 
hrst pinching, when cuj is decreased, of two Stokes lines or turning points which were on 
opposite sides of the real axis for large uji. 

The global frequencies are discretized at order £ and satisfy the interesting property 
of being only (explicitly) dependent on the local dispersion relation (^ in the region 
where the pinching has occurred. Approximations for these modes can easily be obtained 
in the “pinched region” and in the domains of uniform validity of the WKBJ approxi¬ 
mations (white domains of Figures |^), but one must emphasize that approximations are 
not necessarily known or easily accessible everywhere on the real axis. 

^One can point out that nothing can in general be said concerning more stable global modes. 
Furthermore, non-generic cases corresponding to a first pinching of more than two Stokes lines or 
turning points may also give rise to a global mode which is not described by the present analysis. 






The extension in the weakly nonlinear regime has been studied in Chomaz et al. 
(1990), Le Dizes et al. (1993) and Le Dizes (1994) by adding a cubic nonlinearity to 
the right-hand side of equation (|^). In the simplest case where the single destabilized 
linear global mode and its adjoint mode admit uniform WKBJ approximations on the 
real axis, the weakly nonlinear analysis “a la Stuart-Landau” has given the surprising 
conclusion that a physically acceptable solution is only obtained under very restrictive 
conditions (Le Dizes 1994). It is then very likely that there only exists a small class of 
linear global modes which can effectively saturate through a weakly nonlinear process 
when destabilized. This class remains to be determined in order to fully justify the 
application to a physical problem of the linear frequency selection criterion obtained in 
the present article. 

I thank Patrick Huerre and Jean-Marc Chomaz for their helpful advice and kind 
encouragement. 
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